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Abstract 

The spin polarized charge transport is systematically analyzed as a thermally driven stochastic 
process. The approach is based on Kramers' equation describing the semiclassical motion under 
the inclusion of stochastic and damping forces. Due to the relativistic spin-orbit coupling the 
damping experiences a relativistic correction leading to an additional contribution within the spin 
Hall conductivity. A further contribution to the conductivity is originated from the averaged 
underlying crystal potential, the mean value of which depends significantly on the electric field. 
We derive an exact expression for the electrical conductivity. All corrections are estimated in lowest 
order of a relativistic approach and in the linear response regime. 
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I. INTRODUCTION 



The scattering of an unpolarized electron beam at a likewise unpolarized target leads to 
a separation of electrons with different spins 1 . The reason for that effect is the spin-orbit 
interaction originated by relativistic corrections to the Schrodinger equation of the order 
c -2 . As a consequence an electric current in a semiconductor should be accompanied 
by a spin current perpendicular to the conventional current of the charge carriers 2 .. The 
generation of a transverse spin current due to an external electric field called as spin Hall 
effect 3 has attracted much attention recently particularly for realizing spintronic devises^. 
In a series of experiment s^ 1 ^ the effect had been observed in different materials. Recently 
the spin Hall conductivity was measured at room temperature 8 , an observation which give 
the motivation for the present paper. The transport process is considered as a thermally 
activated process which is additionally driven by the applied electric field. This field 
however reveals a direct coupling to the crystal potential in such a manner that also in 
case of a periodic potential a constant current is maintained. Our approach is based on 
the Fokker-Planck equation or Kramers' equation for the semiclassical motion of the charge 
carriers with damping and coupling to a stochastic source like a heat bath. Thus the charged 
particles can overcome permanently, but with a certain probability, potential barriers. Even 
for low temperature one should take into account the coupling of the homogeneous electric 
field to the potential landscape. Actually as discussed below the averaged force acting on 
the particles consists of several parts which point at different directions and contribute to 
the total conductivity. Recently a very promising theory has been proposed^ by studying 
the semi-classical deterministic equation of motion. The theoretical approach describing 
the spin Halle effect goes back to 2 where a mesoscopic equation for the spin density had 
been proposed. A more phenomenological approach based on an anomalous scattering 
mechanism in the absence of spin-flip scattering discussed in 3 . As a consequence anomalous 
currents and a finite spin-diffusion length had been studied in^ where the calculations are 
based on the Boltzmann equation. In a series of papers^^ 3 -"^ a quantum approach for 
the spin Hall conductivity had been discussed. In^ a classical theory is proposed. 
In the present paper we consider an extension of the theoretical approach proposed in&i& by 
including systematically a coupling to stochastic forces and friction terms. Because we are 
in the linear response regime the canonical friction and the strength of the stochastic force 
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are related by the fluctuation-dissipation theorem. The steady state solution of Kramers' 
equation is an appropriate tool to obtain the electrical current. As a novelty we discuss in 
detail the influence of the underlying crystal potential. Unlike^ we get corrections to the 
spin Hall conductivity by the averaged crystal potential. Because the mean value over the 
crystal potential is performed with the stationary solution of Kramers' equation and this 
solution depends on the electric field there occur additional corrections to the conductivity 
in both the linear response regime with respect to the field E and in lowest order c~ 2 of 
an relativistic approach. Whereas the relation for the conductivity under inclusion of the 
crystal potential is an exact result, the real contribution, originated from the potentials, 
has to be estimated by model potentials. 



II. MODEL 

In lowest order of relativistic effects the Hamiltonian has to be supplemented by the 
spin-orbit interaction 

W=^- + [/(r) + Kp-[(TxVU] with k = — %- (1) 
2m 4m 2 c 2 

The quantity a are Pauli operators which are assumed to be constant. Due to the spin-orbit 
interaction the total potential U(r), related to all forces acting on the charged particles, 
appears twofold. In the present paper we incorporate in the potential the intrinsic crystal 
potential U c (r) and the energy of the homogeneous electric field U(r) = U c (r) — eE ■ r. 
Extending the deterministic approach proposed in^ let us include stochastic forces T}{r,t) 
due to other degrees of freedom as imperfections. Additionally, the charged particles are 
subjected to an arbitrary damping force T(r,p). Under that conditions the semiclassical 
equations of motion read 

r = — + kct x W; 

m 

p = - W - kV [p ■ (cr x W)] - r(r, p) + n(t) . (2) 

Because we are interested in the linear response regime the system is assumed to be nearby 
the equilibrium. Therefore it seems to be adequate to assume the stochastic force r) dis- 
tributed by a Gaussian white noise with an arbitrary noise strength D: (rj a (t) r}p{t')) = 
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2 D a p5{ t — t') . The corresponding Fokker-Planck or Kramers' equation for the probability 
density P(r, p, t) 
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(3) 



It is easy shown that the last equation satisfies the principle of detailed balance and exhibits 
an equilibrium solution of the form 

H , 



P e (r,p) ~ exp( 



k R T' 



(4) 



provided the arbitrary damping function T fulfills the relation 

D OH 



k B T dp 



(5) 



From here we conclude immediately 



7II, with II = [p + mK(cr x W(r))] 



(6) 



The quantity II is the canonical momentum whereas the parameter 7 is the damping 
constant or the inverse relaxation time, which appears in the conventional Einstein relation 
D = rwyksT . Due to the linear momentum part in the Hamiltonian ([T]) the 'canonical' 
friction includes likewise a relativistic correction, which leads to an additional term in the 
Drude conductivity discussed below. 



III. STATIONARY CURRENT 

Owing to the electric field energy incorporated into the Hamiltonian Eq. ([TJ, the equi- 
librium distribution Eq. (j2J) may not be appropriate for analyzing the electric current. All 
charged particles would shift to —00 to establish an equilibrium state. Such a state is not 
realized because the equilibrium distribution is not normalizable in that case. Instead of 
the equilibrium solution one needs a steady state solution. Due to the conservation of the 
probability Eq. ([3]) can be written in form of a continuity equation where the probability 
current is a function of both, the spatial coordinate and the momentum j p (p, r, t). Whereas 
the equilibrium distribution P e satisfies j p = a steady state solution P s obeys the weaker 
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condition V • j p = 0. The complete evolution equation for the problem reads now 

dP d d d 2 

1* =-*< AP) -aF (BP > + i V P «" th 

A = — + no- xVU 

m 

B = - W - rcV [p • (a- x VC/)] - T(r, p) (7) 

From here we conclude 

^n) = (B)+m K (A.V(<TX W(r)), ^(r) = (A) . (8) 

To maintain a stationary electric current j the averaged velocity of the charge carriers should 
be fixed. In according to Eq. (JHJ) this requirement is fulfilled by 

j = ne(A>, !<n) = 0. (9) 

The quantity n is the density of the charge carriers. The bracket means now the average 
with the steady state solution P s (r, p; E, k) of Eq. (J7j). It is important to note that all 
averages performed with P s depends on the electric field E, too. Because we are looking 
for a solution linear in the relativistic factor k we consider firstly the non-relativistic case 
k = 0. From Eq. ([8]) combined with Eq. (Q we get 

7<p) = eE-<VC/ c ). (10) 

Since in the non-relativistic = the system offers only one preferential direction, 

namely that one given by the external field E, and furthermore the distribution function of 
the steady state depends on the electric field, the only possibility is that the averaged force 
in Eq. ([TO]) point at the field direction. Hence we set 

(VU e ) =aieE, 

where the dimensionless factor a\ characterizes the influence of the crystal potential. The 
factor ai depends on the concrete form of the potential and will be discussed below for a 
Coulomb potential. Inserting the last relation in Eq. (Q it results in the non-relativistic 

case 

jo = — (l-ax)E. (11) 

7717 

The first part is nothing else as the conventional Drude conductivity with the simplest 
form of a scattering mechanism included in the damping factor 7. The charge carriers will 



follow that law if the influence of the crystal potential is negligible, a situation which should 
be realized at sufficient high temperatures. In case the lattice potential becomes relevant 
an additional scattering mechanism is established and the mobility of a charge particle is 
further reduced leading to a decrease of the conductivity. Now let us take into account the 
relativistic effects in lowest order in k. From Eqs. ([H]),© 
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E--(Wc) 
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where the parameter C is defined by 

d 2 U, 



2nne 2 C (1 — q;i ) . , . 

+ 4 -o-xE, (12) 



C^ v . (13) 



Notice that the factor C differs by a factor e from that defined in^. The last relation 
is obviously valid for an isotropic crystal. As observed in Eq. (|T2l the relativistic effect 
proportional to k yields contributions to the current j pointing toward the field direction E 
and perpendicular to that one in the er x E direction. Hence there exists a third direction 
oriented toward er x (er x E). Because the average over the crystal potential has to be 
fulfilled with the stationary solution of the Kramers equation and this distribution depends 
on the electric field, the most general expression for the averaged force linear in the field E 
reads 

(W c ) = «i eE + Ke[a 2 (c x E) + « 3 (er x (er x E) ) ] . (14) 

Here the parameters a^i =1,2,3 reflect the influence of the crystal potential. They will 
be estimated later. As the final result the electric current consists of three parts 

ne 2 \ ^ ne 2 K«3 , 

j = (1 - ai)E + (er x E) x er] 

2 

+ V±JL[2C{l-oti) -a 2l ]{a x E) . (15) 
In the subsequent section we discuss the conductivity. 



IV. SPIN HALL CONDUCTIVITY 

From Eq. (jl~5l) one obtains the charge conductivity and the Hall conductivity 

2 2 

<Jc= (1-ai , o s = 2Cl- ai -a 27 . 16 
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The parameter C is defined in Eq. ( TT3|) whereas the coefficients are introduced in Eq. ({T41) . 
As already observed in^ the ratio 

_ h 2C(l-ai)-a 2 7 
a c 4m 2 c 2 (1 - ax) 7 

is independent on the concentration of the charge carriers and is determined by the periodic 
potential. The result is in accordance with microscopic models^. The expression for the 
electrical current is different from that obtained previously 9 -. However, if all parameters 
ai are neglected the ratio in Eq. (IT?]) is in accordance with 9 . Apart from the thermal 
activation process as driving force the conventional Drude current is already supplemented 
by a correction term due to the crystal potential, see jo in Eq. (fTTj) . Otherwise the spin 
Hall current is likewise corrected twice by the relativistic damping force T and through the 
influence of the lattice potential. To discuss the differences in more detail let us estimate the 
correction coefficients C, a% . The parameter C is originated by the distribution of charge 
carriers in a cubic lattice. Following the line proposed in^ one can estimate the coefficient 
C using the Poisson equation 

V 2 [7" c (r) = -4vrpe 
where p is the local charge density of the ions. This gives 

C = —Zen (18) 

Here — Ze is the charge and no is the concentration of the ions. Let us notice that the result 
is correct in linear order of the external field. Using the approach proposed below one can 
show that corrections to the parameter C due to the field E occur at first in second order. 
The remaining quantities denoted by a{, may be figured out in a more laborious analysis. In 
finding out the parameters we have to calculate (VU C ) according to Eq. ( 1141 . The problem 
one is confronted with is the fact that the electric field affects the intrinsic crystal potential 
U c (r) in such a manner that a constant current is maintained. In particular, the coupling 
between the external field and the crystal potential becomes important for low excited 
charge carriers which are subjected strongly to the field. In our approach the influence of 
the homogeneous field E on the potential is manifested by the explicit field dependence of 
the steady state solution P s (r, p; E, «). As a direct consequence the averaged force carried 
on the charge particles via the crystal potential, depends on E according to Eq. (fT4"j) . Owing 
to the pure spatial dependence of the crystal potential the reduced distribution function, 

7 



defined by p(r; E, k) — f dpP s (r, p; E, k) seems to be a more appropriate quantity as P s . 
However the basis Kramers' equation ([7j) allows no exact steady state solution P s . Hence 
one is relied on a reasonable approximation for the stationary distribution function. We are 
looking for a solution of Eq. (j7|) in the linear response regime of the external field and in 
lowest order of the relativistic factor k. Let us make the ansatz for the steady state solution 
of Eq. (J7J) in the form 

P 9 (r, p; E, K ) = P (r, p; E) [1 + np ■ f (r, E )] , (19) 

where f = —jKpcr x VU C in lowest order. Inserting this ansatz in Eq. (171) , performing the 
momentum integration and collecting all terms in linear order in E and k one obtains 



(W c ) = Jdr [p(r, E, k = 0) + «p(r, E) • /(r, E) 



W c . 



Here the quantity p(r, E) = J <ippPo( r > P! E) is the averaged momentum proportional to 
the current at the coordinate r for k = 0. Assuming that the reduced distribution function 
p obeys the Smoluchowski equation^ and making the approximation 

«/. *v. \ /|^\ (20) 



, 9x a dxpdx-y I \dx 
we find 

(VC/c) = aieE + ^ ^(cr x E) . 
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Comparing the result with the ansatz made in Eq. (I14p we conclude 

(l-ax)C , . 

a 2 = a 3 = 21 

7 

In lowest oder of a relativistic approach there is no current in spin direction. Inserting this 
result in Eq. ffTTj) one obtains that the ratio of the spin Hall conductivity and the charge 
conductivity is independent on the lattice potential 

a « foe 

(22) 



a c 4m 2 c 2 7 

This expression differs from that obtained in^. by a factor 2. The reason for the discrepancy 
is that in^ the field dependence of the averaged crystal potential is neglected. The spin 
and the charge current is subjected to the same scattering mechanism by the underlying 
isotropic crystal potential. If one includes higher order corrections appearing in Eq. ( 1201 
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then the ratio should be dependent on the potential. Making a similar approach one can 
derive the remaining factor a\. We find in d dimensions 

ai = yd(Lr) 2 / dr( ^ c(r) ~^ 2 ' ^ = vJ drUc{r) (23) 

Notice that the result is valid in first order of k and in the linear response regime. So the 
problem in finding out the influence of the crystal potential is reduced to calculate the mean 
square displacement according to Eq. (123]) . To that aim we have considered a Coulomb 
potential U c = Ze 2 £\ . I r i ~ r j \~ ■ A straightforward calculation yields 

A(Ze 2 ) 2 na _ ( 2e A 



3(k B T)H V V^nk B T , 
where we have assumed a simple cubic lattice with the lattice spacing a. There occurs a 
characteristic energy €a = Ze 2 y/na. For high temperatures €a *C k B T the crystal potential 
becomes irrelevant, i. e the motion of the charge carriers is independent on the underlying 
potential. Otherwise our approach is a classical one which should be restricted to < k B T. 
If the density n of the charge carriers is of the order n ~ a~ 3 then €a — — ■ For a typical 
electronic system the thermal energy is of the order of the Fermi energy k B T ~ ep and 
therefore 

2 

ne 

a c ~ 1 



2e, 



37T 



Even in case that €a — £f there is a significant contribution due to the potential. In 
general the characteristic energy should be smaller than the Fermi energy. Instead of using 
a Coulomb potential one can also consider a cos-potential U c (t) = Uq X^=i(1 — cos^Xj)) . 
For this potential one can calculate the conductivity numerically. The ratio of the 
charge and the spin Hall conductivity versus the rescaled temperature are shown in 
Fig. [H For high temperature the potential becomes irrelevant and the classical Drude 
conductivity is achieved, compare the dashed line in Fig. [H The reason for the weak devia- 
tions of <y s /a c from Eq. (|22p is originated exclusively by the approximation made in Eq. (|20p . 



V. CONCLUSION 



We have presented a thermally driven Drude model under inclusion of the relativistic 
spin-orbit coupling. The charge carriers are driven by a homogeneous electric field and 
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a Gaussian white noise. Additionally the charged particles are subjected to an intrinsic 
periodic crystal potential. Our approach differs from the semiclassical description proposed 
m£ in two respects. Firstly, we start consequently from the stochastic Kramers' equation 
including canonical friction and white noise terms, and secondly we take into account that 
the correction terms due to the averaged crystal potential exhibit a pronounced dependence 
on the the electric field and its orientation with respect to the spins. Our approach 
shows that Spin Hall Effect should also observed for finite temperature as it has reported 
recently^. For high temperatures the particles are not influenced by the potential and the 
conventional Drude conductivity appears. When the temperature is lowered the mobility is 
reduced owing to the potential. In our approach the interplay between the electric field and 
the lattice potential is manifested in the field dependence of the averaged force imposed on 
charged particles. 
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FIG. 1: The charge conductivity in units of the classical Drude conductivity versus the rescaled 
temperature (dashed line) and the ratio a s /a c in units of fiC /Am 2 (?^ (full line) 
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